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Abstract. We perform fully relativistic hydrodynamic simulations of the deceleration 
and lateral expansion of a relativistic jet as it expands into an ambient medium. The 
hydrodynamic calculations use a 2D adaptive mesh refinement (AMR) code, which 
provides adequate resolution of the thin shell of matter behind the shock. We find that 
the sideways propagation is different than predicted by simple analytic models. The 
physical conditions at the sides of the jet are found to be significantly different than at 
the front of the jet, and most of the emission occurs within the initial opening angle of 
the jet. The light curves, as seen by observers at different viewing angles with respect to 
the jet axis, are then calculated assuming synchrotron emission. For an observer along 
the jet axis, we find a sharp achromatic 'jet break' in the light curve at frequencies above 
the typical synchrotron frequency, at tj e t ~ 5.8(i?52/ni) 1 / 3 (6o/0.2) 8 ' /3 days, while the 
temporal decay index a (F v oc t a ) after the break is steeper than — p (a — —2.85 for 
p — 2.5). At larger viewing angles tjet increases and the jet break becomes smoother. 

1 The Hydrodynamics 

The hydrodynamic simulations are performed using a fully relativistic 2D AMR 
code. The jet propagates into a cold and homogeneous ambient medium of num- 
ber density n\ cm -3 . We use n\ = 1, typical of the ISM. We assume an adiabatic 
evolution (i.e. that radiation losses do not influence the dynamics). For the initial 
conditions, we use a wedge of opening angle 9q = 0.2 taken out of the Blandford 
McKee |l| self similar spherical solution, with an isotropic equivalent energy of 
-E52IO 52 ergs, where we use E52 = 1- Since the lateral expansion is expected to be- 
come significant only when the Lorcntz factor of the flow drops to 7 ~ l/0o = 5, 
the initial Lorentz factor of the shock was chosen to be r = \/2j = 23.7. 

Figure ([!]) shows a 3D view of the jet at the last time step of the simulation, 
with color-maps of the proper number density, n' , and proper emissivity. Primed 
quantities arc measured in the local rest frame of the fluid. While the number 
density does not change significantly between the front and the sides of the jet, 
the emissivity is significant only at the front of the jet, and drops sharply at an- 
gles larger than the initial opening angle, 9q. The overall egg-shaped structure is 
very different from the quasi-spherical structure assumed in ID analytic models. 



2 The Light Curves 

The synchrotron spectral emissivity, P' v , , is calculated using the physical condi- 
tions determined by the hydrodynamic simulation. The electrons and the mag- 
netic field are assumed to hold fractions e e and e^, respectively, of the internal 
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energy density, e', while the electrons posses a power law energy distribution, 
N(j e ) oc 7 ( T P for 7 e > j m = [(p- 2)/(p- 1)] (e e e' /n'm e c 2 ). For simplicity, we 
ignore the effects of electron cooling and self absorption, so that P' v , oc j/ 1 / 3 for 
j/ < v' m and P^/ cx iA 1- *')/ 2 for v' > v' m , where u' m is the local synchrotron 
frequency of an electron with j e — j m ; F v is calculated using the formalism of 
, summing over the contributions from the finite 4- volume of the simulation. 
Figure shows the radio light curves seen by an observer along the jet axis 
{&obs = 0). For simplicity, cosmological corrections are not included. The insert 
shows an optical light curve as seen by observers at three different viewing angles 
with respect to the jet axis: 9 b s /9o = 0,1,2. We obtain an achromatic 'jet break' 
in the light curve for v > v m (tj e t) (as predicted by simple semi- analytic models 
|§) at t jet (6 0bs ), where 0mt jet (9 ) « t jet (0) « 5.8(S 52 /n 1 ) 1 /3(0 o /o.2) 8 / 3 days. 
Defining a, (5 by F v oc v^t a , the shape of the break may be approximated by 

F v = F z/ [(t/t jet )- sai + {t/t jet )- sa -\ ~ 1/s , (1) 

where ot\ — a(t <g; tj et ), 012 — ct(t 3> tjet) and the parameter s(9 b s ) determines 
the sharpness of the break, and ranges between s(9o) « 1 to s(0) « 4.5, indicating 
that the break is sharper at smaller o b s - For v < v m (tjet), there is only a 
moderate and more gradual change in a, until the time when u m = v. For 
v > v m we find that 012 is slightly smaller than — p (the value predicted by most 
simple models) for o b s = {a.2 = —2.85 for p = 2.5). 

3 Discussion 

We find that the physical conditions at the sides of the jet are significantly 
different than at the front of the jet, and most of the radiation is emitted 
within the initial opening angle of the jet [9 < 9q; see Figure @]. Therefore, 
the frequently used assumption of a homogeneous jet seems inadequate. For 
v > v m {tjet) we find a sharp achromatic break in the light curve at tj et — 
5.8(£ , 52/ni) 1/ ' 3 (0o/O.2) 8 / 3 days, for 9 a b s — 0, while at larger viewing angles tj et 
increases and the break becomes smoother. For v < v m (tj e t) the change in the 
temporal index a near tj et is more moderate and gradual. The value of a for 
t > tj et , v > v m , 9 bs — is slightly smaller than — p (a = —2.85 for p = 2.5). 
Finally, we note that tj et (0 o bs = 9q) ~ 1.5tj et {9 b s = 0). Since, in order to detect 
a burst in 7-rays we require that 9 b s <> 9q , this may induce an uncertainty of up 
to 15% when deducing the value of 9q from tj e t, unless o b s is well constrained. 

References 

1. R.D. Blandford, C. F. McKee: Phys. of fluids, 19, 1130 (1976) 

2. J. Granot, T. Piran, R. Sari: ApJ, 513, 679 (1999) 

3. J.E. Rhoads: ApJ, 525, 737 (1999) 

4. R. Sari, T. Piran, T. Halpern: ApJ, 519, L17 (1999) 



Light Curves from an Expanding Relativistic Jet 3 




Fig. 1. A 3D view of the jet at the last time step of the simulation. The outer surface 
represents the shock front while the two inner faces show the proper number density 
(lower face) and proper emissivity (upper face) in a logarithmic color scale 



> 
LL 





1 1 


100 GHz 


I 1 

8.46 GHz 


1.43 GHz 




' 

I 

I 

T =i/e r 




10" 2 




^obs=0 


\ \ -2.85 
\ \ l 








\ V 


> 

LL 




e obs= e o r\ 




10" 3 




V 

e obs=2e ... \\ 

\ V 
\ \ 


\ V 

350GHz\ 


10" 4 




\ x 

7=i/e r=i/e \ 




10° 

I 


tdays 10 1 






10° 


^days 


10 1 



Fig. 2. Radio light curves (flux density, F v , in arbitrary units, as a function of the 
observed time in days) for an observer along the jet axis. We use e e = cb — 0.1, p = 2.5. 
The observed times, for an observer along the jet axis, when the Lorentz factors 7 of 
the shocked fluid (vertical dash-dotted line) or F of the shock (vertical dotted line) drop 
to I/80, for an extrapolated spherical evolution, are indicated, (insert) Optical light 
curves for observers at viewing angles 8 b s /8o = 0, 1, 2 with respect to the jet axis. 



